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HEAT TRANSFER TO BODIES IN A HIGH-SPEED RAREFIED-GAS STRAM 


By Jackson R. Stalder, Glen Goodwin, and Marcus O. Creager * 


SUMMARY 

The eiiuilibrium temperature and heat-transfer coefficients for 
transverse cylinders in a high-speed stream oj rarefied gas were 
measured over a range of Knudsen numbers {ratio of molecular- 
mean-free path to cylinder diameter) from 0.025 to 11.8 and 
for Mach numbers from 2.0 to 3.3. The range of free-stream 
Reynolds numbers was from. 0.28 to 203. 

The models tested were 0.0010-^ 0.0050-y 0.030-, 0.05 1-, 
0.080-, and 0.126-inch-diameter cylinders held normal to the 
stream . 

For the case of high-speed flow of a rarefied gas about cylinders y 
the data indicate that fully developed free-rnolecule flow first 
occurs at Knudsen numbers of a 2 )proximately 2.0 and the 
temperature-recovery factor depends primarily upon the Knudsen 
number. For Knudsen numbers greater than 0.2, the tempera- 
ture-recovery factor exceeds a value of unity even though free- 
molecule flow may not be fully developed. 

Over the range of conditions covered in the jyresent tests, the 
Nusselt number is a function only of the Reynolds number if 
the viscosity and thermal conductivity are based upon stagnation 
temperature and the density based upon free-stream conditions. 
For free-molecule flow, the heat-transfer data are well correlated 
by the theory with the accommodation coefficient equal to 0.90. 

INTRODUCTION 

The subject of lieat transfer to or from bodies in a lii^li- 
speed rarefied gas stream has become of importance in 
aeronautics. The recent practical realization of missile flight 
at very high altitudes has stimulated interest in the subject. 
Also, in certain aspects of wind-tunnel testing, notably turbu- 
lence studies in supersonic wind tunnels by means of hot-wire 
anemometers, the eflects of slip may not be negligible. (See, 
e. g., reference 1.) 

The various phenomena encountered in high-speed flow of 
gases are often described in terms of the Reynolds number 
and the Alacli number of the flow. The former parameter 
may be considered to be a measure of the effect of viscosity 
and the latter a measure of the effect of compressibility^ on 
the flow field. 

When considering the flow of rarefied gases, however, an 
additional parameter, the Knudsen number, becomes impor- 
tant. The Knudsen number is defined as the ratio of the 
mean-free molecular path to a characteristic body dimension. 
This parameter, which is a measure of the degree of gas 
rai^faction, may also be construed as defining the importance 

1 Supersedes NACA TN 2438, “Heat Transfer to Bodies in a High-Speed Rarefied-Gas 
Stream” by Jackson R. Stalder, Glen Goodwin, and Marcus O. Creager, 1951. 
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of the effect of molecular motions on the flow field. When 
tlie Knudsen number is small (less than, say, O.OQl), the 
effect of the molecular motions on the flow is negligible and, 
in this regime, the gas can be treated as a continuous medium. 
On the other hand, for large values of the Knudsen number 
(of the order of 10), the effect of molecular motions is all- 
important and the phenomena which occur can be completely 
descrilied in terms of tlie motions of individual molecules — 
in a statistical sense, of course. By virtue of the well- 
established laws of molecular motions, this so-called “free- 
molecule” regime is readily susceptible to analysis and has 
been thoi*oughly examined (references 2 through 8) with 
respect to both aerodynamic and thermodynamic phenomena. 
Experimental work in this j*egime with respect to extez-nal 
flows (as opposed to flow in tubes) is, however, (piite nieagei*. 
The pioneer work of Epstein in 1924 (reference 3) was 
concerned with the measurement of the drag of spheres; the 
tests were confined to very low-speed motions. Recently, 
the drag and equilibrium temperature of a cylinder in 
high-speed free-molecule flow wej*e reported in reference 8. 

The regime of intermediate Knudsen numbers (0.001 
<10) is difficult to deal with analytical!}" and, at this time, 
little success has been attained in attempts to analyze this 
portion of the field of fluid mechanics. Some progress has 
been made, however, by Chapman and Cowling (reference 
9), Schamberg (reference 10), and Grad (reference 11). 

The two general types of approach to the problem used 
by the above authors have been described by Uhlenbeck 
(reference 12). Briefly, the first approach consists of apply- 
ing correction terms to the Navier-Stokes continuum equa- 
tions and boundary conditions, the first approximation 
resulting in a set of higher-order equations — the Burnett 
equations. The other approach is to attempt to solve the 
Boltzmann equation for the molecular velocity distribution, 
which, in effect, results in obtaining correction terms to the 
Maxwell velocity distribution equation. Both of these 
methods are essentially perturbation methods. In the first 
approach solutions are obtained in terms of power series of 
the Knudsen number, while in the second method solutions 
are obtained in terms of the reciprocal of the Knudsen 
number. Both methods will fail for Knudsen numbers of 
the order of unity. 

As in the case of free-molecule flows, little experimental 
work has been done in the intermediate regime for the case 
of high-speed external flows. However, Kane (reference 
13) has reported results of drag tests on spheres in this 
regime. 
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The general purpose of the research described in this paper 
was to study the heat-transfer processes that occur between 
a body and a high-speed rarefied gas stream. More par- 
ticularly, it was hoped that the experimental data would 
yield two important results: first, the proper grouping of 
variables to obtain correlation of heat-transfer data in the 
analytically intractable ‘‘slip-flow’’ regime intermediate 
between the continuum and free-molecule regimes; and 
second, to define approximate values of the parameters 
wdiich delineate the three regimes for the particular con- 
figuration used in these tests. 

Experimental data were obtained on the heat transfer 
from circular cylinders in the three regimes of free-molecule, 
slip, and continuum flows. The data are compared at small 
values of Knudsen number with the data of reference 1 and 
for large values of Knudsen number, with the free-molecule 
analysis of reference 8. 

NOTATION 


Cp 


d 

m 


y(«) 

h 

/„ 

/, 

K 

k 

I 

M 


Nu 

V 

/V 

0 


r 


R 

Re 

s 

T 

n 

T, 

U 

Va 

'^m 


specific heat at constant pressure, foot-pounds per 
slug, °F 

cylinder diameter, feet 
dimensionless function of delined by 

f{s) = Zi(s’^+3) + Z2 

dimensionless function of s defined by (j{s) = 3{Zi~]rZ2) 
gravitational acceleration, 32.2 feet ])er second squared 
heat-transfer coefficient, foot-pounds per scjuare foot, 
second, °F 

modified Bessel function of first kind and zero order 
modified Bessel function of first kind and first order 


Knudsen number 


a>‘'> 


dimensionless 


Boltzmann constant, 5.66X10“^'* foot-pounds per °F 
per molecule 

mean-free molecular path, feet 


Mach number 




I; dimensionless 


Xusselt number 


(t> 


dimensionless 


static pressure, pounds per scpiare foot 


Brand tl number 




; dimensionless 


heat transferred from body, foot-pounds per square foot, 
second 

temperature-recovery factor, defined by 




dimensionless 


gas constant, foot-pounds per pound, °F absolute 
Reynolds number dimensionless 


ss 


molecular speed ratio dimeiisionle 

free-stream static temperature, °F absolute 
temperature of the cylinder, °F absolute 
equilibrium temperature, °F absolute 
stream mass velocity, feet per second 
acoustic speed (-yJygRl"), feet per second 
most probable molecular speed {-^2gRT), feet per 
second 


mean molecular speed 



, feet per second 


Zi dimensionless function of s defined by 


Zj dimensionless function of s defined by 

a accommodation coefficient, dimensionless 

7 ratio of specific heats, dimensionless 
K thermal conductivity, foot-pounds per second, square 
foot, °F per foot 

fx viscosity, pound-seconds per square foot 

V kinematic viscosity square feet per second 

p gas density, slugs per cubic foot 

SUBSCRIPT 

10 refers to the fact that the viscosity and thermal con- 
ductivity were evaluated at stagnation temperature 
while the density was evaluated at free-stream condi- 
tions 

ANALYSIS 

SOME GENERAL RELATIONS BETWEEN VARIABLES 

Before proceeding wdth the analysis, it is desirable to set 
dowTi a few relations which are used tlu*oughout. The 
Knudsen number, = can be related to the Mach and 
Reynolds numbers as follows: 


M_ Vg 

Re U dp 

M 


( 1 ) 


From reference 14, page 147, the viscosity of a rarefied gas 
c-omposed of hard elastic spheres and having a Maxwellian 
velocity distribution can be expressed in terms of the density, 
mean molecular velocity, and mean-free molecular path, 
fx=0A99pvL Thus, very closely. 


Re~ 



( 2 ) 


The molecular speed ratio s, which is a primary variable in 
the analysis of free-molecule flows, is defined as s=UjVm- 
Since v,n=^J'^{jRT and Va=^y(jR7\ we have a relation 
between .s and \I 

(3) 

HEAT TRANSFER FROM CYLINDERS IN FREE-MOLECULE FLOW 

Heat transfer from a circular cylinder oriented trans- 
versely to a gas stream has been treated, for the case of free- 
molecule flow, hi reference 8. It is useful, however, to 
examine the final equations derived therein wdth a view to 
putting them in such form as will allow plotting of the 
experimental data on a basis comparable with continuum 
data. 


HEAT TRANSFER TO BODIES IN A HIGH-SPEED RAREFIED-GAS STREAM 


3 


The equation for heat transfer from a cylinder to a diatomic 
gas stream is shown in reference 8 to be (neglecting radiation 
terms) 


Tc 

T 


gis)-f(s) 


pVmOC 


Q=0 


(4) 


One can define a heat-transfer coefficient per unit area as 



Te 


(5) 


where Te is the equilibrium temperature assumed l)V the 
cylinder in the absence of heat transfer. The quantity Te 
can be obtained from equation (4) by setting the term Q 
equal to zero in which case Tc equals Te. Thus, 


rp rp 

g{s) 


( 6 ) 


Then, substitution of equations (4) and (6) into equation 
(5) yields 


,_pVmOig(s) 

27t^/2T 


(7) 


If use is made of the relation pjT=gpRj equation (7) can be 
rewritten as 

^'“■"2^1: <*) 

Substituting numerical values for the constants, for the 
case of nitrogen gas, we have, 

Nu = 0.02G4aReFr^ (9) 

The term g(s)ls is nearly constant for values of s greater than 
2 and approaches a value of 6^^ asymptotically. The 
limiting values are 



Then for values of s>2 we have, very nearly, 

NyF=0.2S3anePr ( 10 ) 

For gases, the Prandtl number, Pr, is nearly constant over 
a wide temperature range and, for diatomic gases, has a 



Figure l.— Ames low-density wind tunnel. 
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value of about 0.72 at 0° F. Substitution of this value for 
Pr into equation (10) yields the final equation for heat trans- 
fer from a transverse cylinder in a high-speed (.9>2) free- 
molecule flow field 

Xu=0.2inane ( 11 ) 

DESCRIPTION OF EQUIPMENT AND EXPERIMENTAL 
TECHNIQUE 

THE WIND TUNNEL 

The Ames Laboratory low-density wind tunnel, in which 
the present tests were conducted, is an open-jet, nonreturn- 
type tunnel and is cl escribed in detail in reference 8. A 
sketch of the major components of the wind tunnel is shown 
in figure 1 . 

The nozzles used in the jirc^sent test, however, diflered 
from those descril)ed in reference 8 and foi* that reason will 
be described in detail here. The nozzles used were axially 
symmetric, iK inclu's in exit dianietcT, and had porous walls 
to allow part of the boundary layer to be removed by suction. 
The nozzle contours were determined by calculating the 
shape of the inner core (assumed inviscid) by the method of 
characteristics, as outlined in reference 15, and applying a 
boundary-laycu* (‘orrection to the contour. The boundary- 
layer correction was calculated using the von Karnuin 
moment um-int(‘gral tlu'ory given by Schaaf in reference 16. 
The correct icm actually auplical to the nozzle contours 
corresponded to that calculated for the highest ])ressure level 
used. 

The nozzles were constructed by stacking 2-inch-diameter, 
thin metal shims to the desired length. The stack of shims 
was held together by four through-bolts and was machined 
to the calculated internal contour. The nozzle walls were 
made porous by discarding every other solid shim and re- 
placing it with small spacers around each of the four through- 
bolts. The shim thickness was varied along the length of 
the nozzle in order to olitain a reasonably uniform suction 
velocity normal to the wall. The nozzle boundary-layer 
removal was accoui])lished by connecting a chamber around 
the nozzle to one of the main vacuum pumps. The porous 
nozzle is shown diagrammatically in figure 2. 



Figure 2. — Cutaway view of porous nozzle. 


The use of boundary-layer removal resulted in a nozzle 
which produced an essentially constant Mach number ov^er 
the center third of its diameter in contrast with the solid 
nozzles described in reference 8 which produced essentially 
a parabolic Mach number distribution. Also the use of 
boundary-layer removal permitted the pressure level of the 
wind tunnel to be changed by a factor of 2}i, while the Mach 
number varied only by about ±2 percent. A Uqiical Mach 
number distribution obtained with and without boundarv- 
la^^er suction is shown in figure 3. 



Figure 3.— Variation of Mach number with distance from center of stream with and without 
boundary-layer suction. Stream static pressure 115 microns of mercury absolute. 


The three nozzles used in the tests produced Mach numbers 
of 2.00±1K percent, 2.50±2 piwcent, and 3.15±2 percent 
over a range of pressure levels from 80 to 200 microns of 
mercury absolute in the center lialf-inch of the jet. 

MODELS TESTED 

ddie models tested were right circular cylinders held 
normal to the air stream by co])per rods fastened to the ends 
of the cylinders. Six models were used in the tests; their 
diameters, test lengths, and types of construction are given 
in the following table: 


Model 

No. 

Material 

Diame- 
ter (in.) 

Test 

length 

(in.) 

Remarks 

1 

Pt-Ni alloy 

0.0010 

0. 434 

Solid. 

2 

Pt 

.0050 

.480 

Solid. 

3 

Fe — Constantan 

.030 

.50 

Butt-welded thermocouple. 

4 

Pt film on Pyre.x glass. . . 

.080 

.530 

Solid. 

5 

Ni 

.051 

.:J80 

fNi wire wound on ceramic 

G 

Ni 

.126 

.3G8 

\ tube. 


Models 1 through 4 had guard heaters wound on the copper 
support rods to allow the temperature of the ends of the 
models to be adjusted in order to prevent end losses. Ther- 
mocouples were soldered to each end support in order that 
the model end temperature could be measured. Figure 4 
shows the arrangement of the model and its supports in the 
stream. Models 5 and 6 were constructed by wrapping 
0.003-inch-diameter insulated nickel wire on a ceramic tube. 
Voltage leads were soldered to each end of the test length and 
brought out through the ceramic tube to the ends of the 
model. The model was then dipped in a thin solution of 
enamel to obtain a smooth surface. These models were made 
long enough to completely span the nozzle. 
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Base 


Figure 4— Diagrammatic sketch of model and supports. 


Model 3 was a butt-welded iron-constantaii tliermoeoiiple 
with the junction located in the center of the test length. 

Model 4 was constructed by coating a Pyrex glass rod with 
a very thin layer of platinum. The model was mounted as 
shown in figure 4 by soldering the copper rods directly to the 
platinum coating. 

The temperature of all the models, except number 3, was 
determined from their resistance which was measured by 
comparing the voltage drop across the model to the voltage 
drop across a standard resistor connected in series with the 
model. These voltages were measured with a direct-current 
potentiometer and galvanometer havmg a least count of 1 
microvolt. The temperature-resistance characteristics of 
each model were determined by calibrating the model in a 
constant temperature bath over the range of temperatures 
encountered in the tests. The temperature of model 3 was 
determined by measuring, with a potentiometer, the thermo- 
couple voltage generated at the junction. 

TEST PROCEDURE 

An impact-pressure tube, a static-pressure tube, and the 
test model were mounted on a common support which could 
be moved tlu*ough the test area of the jet. The wind tunnel 
was then started and the impact and static pressures were 
measured over the test area of the nozzle. The static- 
pressure tube has been described in reference 8 and the im- 
pact-pressure tube was identical with the 0.150-inch-diameter 
impact-pressure tube described in reference 17. Pressure- 
measuring techniques are also described in reference 8. 44ie 
impact-pressure corrections due to viscous effects, given in 
reference 17, were used in order to calculate the stream Mach 
number. 

After the stream conditions had been determined, the test 
model was lowered into the stream and the equilibrium tem- 
peratures of the model test length and ends were measured. 
In the case of models 1 and 2, the guard heaters were used to 
raise the model end temperatures until they were equal to the 
indicated test-length temperature. Heat was not supplied to 
the ends of the other models for the measurement of equi- 


librium temperature as, for the most part, the equilibrium 
temperatures were, on the average, 10° F below the stagna- 
tion temperature of the stream. 

After the equilibrium temperature of the model was 
measured, the model was heated by passing a direct current 
from a batteiy through the model. For models 1, 2, and 4, 
the end temperatures were adjusted until there was no 
temperature difference between the ends and the model test 
length. When this point was reached, the voltage across the 
test length of the model and across the standard resistor in 
series with the model was read on the potentiometer. The 
heat sup])lied to the model was then the product of the 
voltage drop across the test length of the model and the 
current through it. 

Model 3, the butt-welded thermocouple, was not used in 
the heat-transfer tests. 

Independent control of the end temperatures was not 
possible with models 5 and 6; however, these models were 
constructed in such a manner that the resistance to beat 
flow along the models was large and the test length was 
confined to the center % inch of the models, thus minimizing 
the effect of end losses. 

The radiation heat loss from the models was determined 
by measuring the heat loss over a range of model tempera- 
tures, with no flow passing through the tunnel, as the test- 
chamber pressure was decreased toward zero. As it was not 
possible to obtain a test-chamber pressure below 1 micron of 
mercury absolute, the actual radiation heat loss was obtained 
by extrapolating the measured heat loss to zero pressure. 
Models 1 and 6 were later placed in a chamber in which the 
pressure could be maintained at 0.1 micron of mercury ab- 
solute and the results checked the extrapolated data. 

The resistance of the models at room temperature was 
checked from day to day and, if a resistance changed enough 
so that the temperature of the model was in doubt by more 
than 2° F, the model was either recalibrated or replaced. 

The tests were conducted at three nominal Mach numbers, 
2.0, 2.5, and 3.15, and over a range of pressures from 80 to 
200 microns of mercury absolute. 

The test gas used was commercially pure dry nitrogen. 

A summary of the test data is shown in table I. 

RESULTS AND DISCUSSION 

TEMPERATURE RECOVERY FOR AN INSULATED BODY 

The usual starting point for a set of heat-transfer experi- 
ments at high airspeeds is the determination of the equilib- 
i*ium temperature of the body; that is, the temperature as- 
sumed by the body in the absence of heat transfer. In the 
case of a flat ])late in continuum flow, several analyses are 
available for each of the two cases of laminar and turbulent 
boundary layers. In these analyses, the equilibrium tem- 
perature is expressed in terms of ^‘recovery factor,’^ r, which 
can be implictly defined as 

A summary of the available experimental and theoretical 
results is given in reference 18. In brief, the theoretical 
results indicate the recovery factor for laminar boundary 
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TABLE I.— TEST DATA 


Model 

No. 

Mach 

No. 

Static 
pressure 
(Ib/sq ft) 

Stagna- 
tion tem- 
perature 
(°F) 

T. (°F) 

Tc (°F) 

> 

1 

Net heat 
transfer 
from model 
(watts X 10”3) 

Model 

No. 

Mach 

No. 

Static 
pressure 
(Ib/sq ft) 

Stagna- 
tion tem- 
perature 
('^F) 

T, (°F) 

Tc (®F) 

AT= Tc-T, 
(°F) 

Net heat 
transfer 
from model 
(watts X 10-3) 

1 

1.96 

0.280 

69.0 

152.5 

303 

150.5 

3. 03 

2 

3. 19 

0.330 

63.0 

100.5 

310.3 

209.8 

25.9 


2.01 

.529 

71.0 


313 

«160.5 

5. 95 


3. 15 

.437 

67.7 

101.3 

308.0 

206.7 

32.7 


1.97 

.377 

68.5 


308 

“ 161. 8 

4. 32 


3. 12 

.547 

68.5 

96.0 

307.5 

211.5 

41.0 


1 95 

331 

68 7 

148 5 





3. 17 

.378 

64.3 

105.0 





1 99 

* 41^ 

71 3 

148.0 














378 

71 6 

300 

« 146. 5 

4. 83 

3 

3. 11 

.274 

70.4 

60.6 






*378 

71 6 


400 

» 246. 5 

8. 10 


3. 15 

.395 

67.4 

56. 5 






422 

71 3 


301. 5 

" 149. 4 

5.63 


3. 04 

.503 

68.2 

57.8 





2. 48 

. 422 

71. 3 


401. 5 

» 249. 4 

9.80 


3. 17 

.546 

68.5 

57.6 





2 49 

. 468 

68. 7 


300 

» 153.3 

6. 20 


2. 46 

.399 

68.5 

59. 2 





2. 49 

. 468 

68. 7 


403 

«256.3 

10.65 


2.62 

.411 

70.5 

60.7 





2 58 

. 503 

71. 0 


:ioi. 5 

“ 152. 4 

6. 73 


2. 40 

.250 

68.0 

60. 4 





2 . 58 

" 503 

71. 0 


400 

»250. 9 

11. 20 


2. 59 

.517 

68.5 

57.3 






2. 43 

. 341 

72. 3 


302 

« 146. 8 

5. 17 










2.43 

‘.341 

72.3 


404 

o 248. 8 

8.56 

4 

1.95 

.346 

^.5 

66.0 

^2 

166 

106.5 


2.50 

.338 

67.5 


302 

» 152. 3 

4.88 


1.98 

.418 

66.0 

66.0 

223.5 

157.5 

126.8 


2.50 

.338 

68.7 


404 

« 252. 9 

8. 00 


2. 11 

.564 

65.7 

65.7 

203.5 

137.8 

132.3 


2. 52 

. 457 

66. 7 


300 

“ 155. 5 

6. 40 










2. 52 

!457 

66! 7 


398 

“253.5 

10. 47 

5 

1.98 

.350 

69.0 

n.5 

200 

128.5 

45.3 


2. 47 

.376 

65. 0 

146. 6 





2.01 

.504 

71.0 

65.5 

200.3 

134.8 

65.2 


2. 59 

. 493 

68. 7 

144.0 





2. 06 

.239 

70.3 

79.3 

197.2 

117.9 

31.8 


2 44 

. 304 

66. 4 

153. 0 













2. 55 

. 437 

70.3 

148. 5 




6 

2.01 

.324 

66.5 

64. 1 

198.4 

134.3 

85.5 


k 2. 25 

! 180 

72.0 

173.5 

301.5 

128 

1.88 j 


2.01 

.267 

68.2 

68.1 

205.4 

137.3 

76.3 


2. 44 

.224 

65.0 

160.5 

300 

139. 5 

2.49 


2. 04 

. 449 

65.0 

58.4 

176.9 

118.5 

97.8 










2. 04 

.285 

68.7 

69. 1 

186. 1 

117.0 

67.0 

2 

2.47 

.343 

70.3 


301.6 

» 168 

17.0 


2. 05 

.507 

64.3 

57.2 

193.6 

136.4 

123.3 


2.60 

.503 

67.5 


305 

“191.1 

27. 1 ! 


3. 18 

.296 

67.0 

59.0 

197.0 

138.0 

123.7 


2. 56 

.492 

71.5 


302.5 

“ 177. 5 

21.8 1 


3. 36 

.332 

65.0 

53.1 

184.5 

131.4 

146.5 


2. 47 

. 355 

74. 0 

132. 7 



1 


^2.80 

.342 

70.0 

56.4 

178. 0 

111.6 

136.7 


2. 39 

. 304 

68. 9 

136. 7 





^2. 65 

.272 

73.0 

57.5 

181.7 

124.2 

108.2 


*>' 

. 388 

69. 4 

129.0 





f>2. 49 

.217 

74.0 

62.0 

176.0 

114.0 

80.7 


z. 

2. 57 

. 458 

69. 7 

116. 0 





6 2. 51 

. 139 

68.5 
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128! 0 
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6 2. 39 

.087 
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71.0 
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39.6 
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.315 
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120.3 

306.5 

186.2 

24. 0 


6 2.53 

.208 

72.4 

59.7 

175 

115.3 

84.0 


“ T, obtained from figure 5. 

*> Data taken in University of California tunnel. 


layers to be independent of both Mach and Reynolds num- 
bers and to have a value of All the turbulent-bound- 

ary-layer analyses assumed incompressible flow, hence Mach 
number effects do not appear. The analyses, in general, in- 
dicate that the recovery factor for turbulent boundary layers 
on a flat plate has a value close to with some analyses 
showing a minor dependence of recovery factor on the 
Reynolds number. 

For the case of a cylinder with axis normal to the stream, 
little theory exists and recourse to exi)eriment is required. 
Several experiments have contributed information concern- 
ing the equilibrium temperature of transverse cylinders in 
high-speed air streams, notably those described in references 
1, 19, and 20. There is considera))le disagreement among 
the data cited and apparently a definitive experiment with 
a systematic variation of Reynolds and Mach number^ has 
yet to be made. 

The equilibrium temperature for the case of cylinders in 
free-molecule flow has been detennined both theoretically 
and experimentally in reference 8. Good agreement was 
obtained between the experimental results and the somewhat 
surprising theoretical prediction that the cylinder equilib- 
rium temperature would exceed the total or stagnation 
temperature of the gas stream. This phenomenon, of 
course, is in direct contrast to the corresponding phenomenon 
which occurs under continuum-flow conditions where an 
insulated body can have, at most, a temperature equal to 
the stream total temperature and normally does not even 
attain this temperature due to outward heat flow in the 
boundary layer. The anomaly can be explained, however, 
by a consideration, in the case of free-molecule flow, of the 
magnitudes of the incident and re-emitted total velocity 
resulting from combination of the stream mass velocity and 
the random thermal velocities. AVhen the total velocity 


term is squared there results a term, which is twice the 
scalar product of the vector mass velocity and the vector 
thermal velocity. This scalar product term effectively 
increases the apparent internal energy of the gas from the 
continuum value of 3/2 kT per molecule to a value varying 
from 2 kT to 5/2 kT depending upon the speed and orientation 
of the body. The apparent internal energy becomes 5/2 kT 
for high speeds and body angles of attack greater than zero. 
This is just equal to the internal plus potential energy per 
molecule of a cube of continuum gas. Therefore, the energy 
incident on the body in free-molecule flow becomes equal to 
that of a continuum for large values of the molecular speed 
ratio for surfaces inclined at angles of attack greater than 
zero. The molecular energy which is re-emitted from the 
surface is assumed to be equal to the energy of a stream 
issuing effusively from a Maxwellian gas, in equilibrium at 
some yet unspecified temperature, into a perfect vacuum. 
For the case of an insulated body, the temperature of this 
gas is that of the body. The energy of the re-emitted stream 
calculated in this manner is equal to 2 kT^^ per molecule. 
The corresponding energy for a continuum gas at the same 
temperature is 5/2 kTe per molecule if both thermal and 
potential energies are considered. For a given re-emitted 
stream temperature, a smaller amount of energy per mole- 
cule is transported from a body for the case of free-molecule 
flow than is transported in the case of continuum flow. It 
is clear that, if the same total amount of energy, namely 
the incident energy, is to be removed in each case, the 
effusive stream temperature and, hence, the body temperature 
for an insulated body must be higher for free-molecule flow 
than for continuum flow. 

The equilibrium temperature for a cylinder in a free- 
molecule flow field has been shown (equation (6)) to be a 
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function only of the molecular speed ratio, 6‘, or the Mach 
number, AL (See equation (3).) 

Because of the (lifferent ])arameters u])on which the re- 
covery factor or equilibrium temjierature depends in the 
several regimes, it is difficult to plot all the equilibrium- 
temperature or recovery-factor data on a single graph. 
However, for the case of free-molecule flow, the theory ])re- 
dicts that the eciuilibrium-temperature ratio is a function of 
the molecular-speed ratio only. Hence, if we should plot 
the ratio of the measured equilibrium-tem])erature ratio 
TejT to the theoretical equilibrium-temi)erature ratio/(.sO/t/('S*) 
as a function of Knudsen number, the deviation of the ratio 
from a value of unity should define the lower limit of free- 
molecule flow in terms of the Knudsen numl)er. These data 
are shown in figure 5 and indicate that the Knudsen numf)er 
must have a value of approximateh^ 2 in order that fully 
developed free-molecule flow exist for a transverse c}dinder. 

The recovery-factor data are shown in figure G, ])lotted as 
a function of Knudsen number. It can be seen from this 
figure that the data, in the range of Knudsen numbers from 
0.02 to 2.0, are correlated by Knudsen mimi)er alone with no 
S 3 "stematic Mach number effects shown. It is further evident 
that, for values of Knudsen number greater than 0.2, the 
recoveiy factor exceeds a value of unit}’, presumably due to 
the development of free-molecule-flow effects. Also, it can 
be seen that for values of Knudsen number less than 0.2, 
the recovery factor has a constant value indicating an inde- 
pendence from Reynolds number effects. (See equation (2).) 


HEAT-TRANSFER-TEST RESULTS 

The jn-esentation of heat-transfer data for these tests is, 
as with the recovery-factor data, beset with difficulty due to 
lack of a common correlating parameter. The results of the 
heat-transfer tests are shown in one form in figure 7 in which 
Nuiq is plotted as a function of Rei^. The subscript 10 refers 
to the fact that the viscosity and thermal conductivity were 
evaluated at tunnel stagnation temperature, while the density 
was evaluated at free-stream conditions. The use of Nuiq 
and /i^io was suggested by Kovasznay and Tormarck, who 
found that this particular choice of conditions for evaluation 
of the gas pro])erties eliminated the need for including the 
Mach number as an additional variable. The data of 
Kovasznay and Tormarck (reference 1 ) are shown as a dotted 
line which was computed from the expression 

0.580 0.795 

which is an empirical equation resulting in the best fit of 
tlieir data. All the heat-transfer data obtained during the 
present tests are correlated by the bi’oken line with an 
average deviation of ±6 percent. This line was faired 
through the data by the method of least s(|uares and may be 
represented by the relation 

iV?/io=0.132 (12) 

Tlie solid line represents equation (11) of this paper with the 
accommodation coefficient a having the value of 0.9. This 
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value of a was chosen to give the best fit of the data obtained 
in the free-molecide range of Knudsen numbers. 

At the other extreme range of Reynolds numbers, how- 
ever, it is seen that the data obtained in the Ames Laboratory 
low-density wind tunnel indicate Nusselt numbers substan- 
tially below those measured by Kovasznay and Tormarck. 
When these data were obtained, it was suspected that the 
discrepancy was due to blockage effects since the model 
(}^-inch diameter) was large compared to the diameter 
{){ inch) of that portion of the stream suitalde for testing. 
Consequently, some of the tests were repeated using the 
l^niversity of California low-density wind tunnel (see refer- 
ence 21 for description) which has a considerably larger 
(approximately 3-inch diameter) testing region, ddie higher 
Reynolds number data obtained in the Ihiicersitv of Cali- 
fornia tunnel again indicated Nusselt numbers below those 
obtained by Kovasznay^ and Tormarck although higher than 
those obtained in the Ames tunnel. The fact that the 
present low-density data were obtained in two different wind 
tunnels and (in tl\e case of tl\e Ames test) with three different 
models, while the data of Kovasznay and Tormarck are 
corroborated by the data of Lowell (reference 20), indicates 
that the difference in Nusselt numbers is real and is not due 
to experimental error. ^ 

The discrepancy may be due to one or a combination of 
several effects. First, there is the probability that a differ- 
ence in tunnel turbulence level existed between the tests of 
Kovdszna}' and Tormarck and Lowell, which were conducted 
at relatively high pressures, and the present tests which 
were conducted at extremely low pressures. It appears 
likely that the turbulence level in a low-density nozzle woidd 
be low due to the large viscous effects present.^ 

It is also probable that a correction for the effect of vis- 
cosity should be applied to the readings of the static pressure 
tube, as was done in the case of the impact pressure tube 
readings. The value of this correction (which would affect 
the calculated Mach and Reynolds numl)ers) is unknown 
at present. 

It is apparent that additional tests are required to resolve 
the discrepancy. 

CONCLUSIONS 

The following conclusions can be drawn from the data for 
transverse cylinders presented in this paper: 

1. Fully developed free-molecule flow occurs for Knudsen 
numbers of approximately 2 and higher. 

2. For Knudsen numbers greater than 0.2, the tempera- 
ture-recovery factor exceeds unity even though free-molecule 
flow is not fully developed. 

3. For the case of high-speed flow of a rarefied gas about 
cylinders, the temperature-recovery factor is primarily 
dependent upon Knudsen number. 

4. Over the range of conditions covered in the present 
tests, the Nusselt number is a function only of the ReAmolds 

2 It can be observed that the heat-transfer data obtained in the University of California 
tunnel with the 0.001-inch-diameter model are in good agreement with the data obtained in 
the Ames tunnel; howev^er, it should be pointed out that a shift in the temperature calibration 
of this model resulted in an uncertainty of the order of =fcl0 percent in the data obtained in 
the University of California tunnel. 

® In passing, it is noted that McAdams (reference 22) reported large changes in heat-transfer 
rates from cylinders due to variation of the turbulence level of the air stream although insuffi- 
cient details of the original e.xperimental work were given to enable any conclusions to be 
drawn that would be pertinent to the present discussion. 


number if the viscosity and thermal conductivity are based 
on stagnation temperature and the density based on free- 
stream conditions. 

5. In free-molecule flow, the heat-transfer data are well 
correlated by the theory with the accommodation coefRcient 
equal to 0.9. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 14 , 1951 . 
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